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TO ALL WHOM IT MAY CONCERN: 

Be it known that we, Robert G. Waggener, of P.O. Box 3960, Highway 35 N, 
Natalia, Texas 78059-2606, and Elizabeth Wieckowska, of 5943 Heather View, San 
Antonio, Texas 78249, each a citizen of the United States of America, have invented 
certain new and useful improvements in: 



THREE COMPONENT X-RAY BONE DENSITOMETRY 



of which the following is the specification. 
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THREE COMPONENT X-RAY BONE DENSITOMETRY 

CROSS-REFERENCE TO RELATED APPLICATION 

5 The benefit of U.S. Provisional Patent Application Serial No. 60/428,300, filed 

November 22, 2002, is hereby claimed, and the specification thereof incorporated 
herein in its entirety by this reference. 

BACKGROUND OF THE INVENTION 

10 "Osteoporosis is a systemic skeletal disease characterized by low bone mass and 

macroarchitectural deterioration of bone tissue, which leads to bone fragility and 
susceptibility of fracture" (Blake, et ai, "The evaluation of osteoporosis: dual energy x- 
ray absorptiometry and ultrasound in clinical practice." Martin Dunitz Ltd 1999.). 
According to the National Osteoporosis Foundation one in two women and one in eight 

15 men over the age of fifty will have an osteoporotic-related fracture in their lifetime. 
Eight million American women and two million men have osteoporosis and millions 
more have low bone density. The use of noninvasive bone densitometry has resulted in 
significant improvement in the early detection of osteoporosis. 

Dual x-ray absorptiometry (DXA) is the most often used method in the 

20 measurement of bone mineral density because of the low cost, low radiation dose and 
the possibility of measurement at the anatomical site where the fractures most 
commonly occur (femur, spine). (See Genant, et aL 9 "Noninvasive assessment of bone 
mineral and structure: state of the art." J Bone Miner Res. 11: 707-30 (1996)). By 
using only two energies, the dual x-ray absorptiometry method is limited to 

25 distinguishing tissue composition to two components (muscle and bone or muscle and 
adipose tissue). To correct for different tissue composition (varying proportions of lean 
and adipose tissue), the fat fraction is calculated over the soft tissue part of the patient 
where no bone is present. The assumption is made that adipose tissue is uniformly 
distributed and the soft tissue composition adjacent to the bone is the same as the soft 
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tissue overlying the bone. However, in fact, fat is not uniformly distributed and 
therefore, this assumption can introduce errors of greater than 10% in a significant 
percentage of the patient population. (See Svendsen, et ah "Impact of soft tissue on in- 
vivo accuracy of bone mineral measurements in the spine, hip and forearm: a human 
5 cadaver study." J Bone Miner Res 10: 868-73, (1995); Tothill, et ah "Errors due to 
non-uniform distribution of fat in dual x-ray absorptiometry of the lumbar spine." Br J 
Radiol 65:807-13, (1992); Farrell, et ah "The error due to fat inhomogenity in lumbar 
spine bone mineral measurements." Clin. Phys. Physiol. Meas. 10 57-64, (1989)). 

It would be desirable to provide a method and apparatus that can, from the 
10 energy distribution of the x-ray spectra, incident and transmitted through the patient, 
determine the areal densities of the three tissue components: bone, fat and soft tissue. 
The present invention addresses this need and others in the following manner. 

SUMMARY 

15 The present invention relates to deriving the in vivo three-component areal 

density composition of a patient along a transmitted x-ray beam by raster-scanning a 
collimated x-ray beam across a bony region of the patient, deriving the x-ray energy 
spectrum at each point, and using matrix equations with a priori spectra information to 
solve for the relative areal density of soft tissue, fat and bone. 

20 It is to be understood that both the foregoing general description and the 

following detailed description are exemplary and explanatory only and are not 
restrictive of the invention, as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 
25 The accompanying drawings illustrate one or more embodiments of the 

invention and, together with the written description, serve to explain the principles of 
the invention. Wherever possible, the same reference numbers are used throughout the 
drawings to refer to the same or like elements of an embodiment, and wherein: 
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Figure 1 illustrates an apparatus for determining areal densities of bone, fat and 
soft tissue in a patient. 

Figure 2 shows some exemplary curves from an extensive database of measured 
spectral curves for in-air, i.e., unattenuated, x-ray beam and an x-ray beam attenuated 
5 by a phantom simulating body composition of bone, muscle and fat. 

Figure 3 shows exemplary curves of the logarithm of the ratio (Io/I) at the 
photon energy range, 36.5 to 56.5 keV derived from x-ray spectra similar to those 
shown in Figure #2 for both calculated and measured x-ray spectra. 

Figure 4 shows 1 1 exemplary x-ray spectra from 40 kVp to 140 kVp with a 0.5 
10 mm aluminum (Al) added filter. These spectra are drawn from the database which was 
constructed from measured transmission curves. 

Figure 5 shows transmission (attenuation) curves of each of the spectra shown 
in Fig. 4 using an added aluminum (Al) filter. 

Figure 6 is an excerpt of computer code illustrating an iterative process that 
15 takes the sum of the individual solutions and estimates the areal density of the 
individual components. 

Figure 7 shows the composition and fraction by weight of the constituent 
elements provided by the manufacturer for the three types of tissue used in the study. 

Figure 8 shows the pulse height distribution of the photons of the primary 
20 spectrum obtained with an acrylic water container and acrylic in the x-ray beam for 87 
kVp, 0.01 mA tube current and 500 second duration. 

Figure 9 shows the pulse height distribution of the photon spectrum for the 
container after filling the container with 1 5 cm of water. 

Figure 10 shows values of areal densities and the compositions of the phantoms. 
25 Figure 11 shows calibration spectra of Co-57 and Am-241 sources used for the 

calibration of the spectroscopy system. 

Figure 12 shows the energies of the Co-57, Ba-133 and Am-241 sources plotted 
versus channel number. 
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Figure 13 shows ten x-ray spectra taken at 85 kVp, 0.01 mA and 50 second 
duration plotted on the same graph. This was done to demonstrate the reproducibility 
of the x-ray spectra measurements. 

Figure 14 shows the mean values of number of photons, standard deviations and 
5 coefficients of variation of the ten x-ray spectra illustrated in Figure #13. 

Figure 15 shows ten composite x-ray spectra curves taken at 85 keV, 0.01 mA 
and 500 second duration plotted on the same graph. 

Figure 16 shows the mean values of number of photons, standard deviations and 
the coefficients of variation of the x-ray spectra shown in Figure #15. 
10 Figure 17 shows theoretically derived values of the mass attenuation 

coefficients for muscle, adipose tissue and bone simulating phantoms plotted as the 
function of energy in the energy range of interest. 

Figure 18 shows the pulse height distribution of primary and transmitted spectra 
through 5 cm of muscle simulating phantom. The data was obtained at 85 kVp, 0.011 
15 mA and 500 second duration. 

Figure 19 shows experimentally obtained mass attenuation coefficients for 
muscle compared with the theoretically determined values calculated from the pulse 
height distribution photon spectra in Figure #18 in the energy range of 36.5-57.5 keV. 

Figure 20 shows a comparison of experimentally and theoretically derived mass 
20 attenuation coefficients of the muscle simulating phantom. 

Figure 21 shows the pulse height distribution of the photons of the primary and 
transmitted spectra through the 3 cm of adipose tissue simulating phantom. The x-ray 
spectra were generated at 85 kVp, 0.01 1 mA and 500 second duration. 

Figure 22 shows experimentally obtained mass attenuation coefficients for the 
25 adipose tissue simulating phantom compared with the theoretically determined values. 

Figure 23 shows a comparison of experimentally and theoretically derived mass 
attenuation coefficients of the adipose tissue simulating phantom. 



197311 

4 



ATTORNEY DOCKET NO.: 21105.0005U2 

Figure 24 shows the pulse height distribution of the photons of the primary and 
transmitted spectra through 2 cm of bone tissue simulating phantom obtained at 86.5 
kVp, 0.01 mA and 500 second duration. 

Figure 25 shows experimentally obtained mass attenuation coefficients of the 
5 bone tissue simulating phantom compared with the theoretically derived values. 

Figure 26 shows a comparison of the experimental and theoretically derived 
mass attenuation coefficients of the bone tissue simulating phantom. 

Figure 27 shows the pulse height distribution of the photon spectrum 
transmitted through 20 cm muscle, 3 cm adiopose and 1 cm bone simulating phantom. 
10 The spectrum was generated at 82 kVp, 0.343 mA and 750 second duration. 

Figure 28 (Phantom 1) shows the comparison of calculated values of the sum of 
the products of the mass attenuation coefficients and areal densities for three 
components of a first phantom and measured fractional logarithmic transmissions. 

Figure 29 shows the pulse height distribution of the photon spectrum 
15 transmitted through 15 cm muscle, 3 cm adipose and 1 cm bone simulating phantom. 
The spectrum was generated at 82 kVp, 0.152 mA and 750 second duration. 

Figure 30 (Phantom 2) shows a comparison of calculated values of the sum of 
products of the mass attenuation coefficients and areal densities and measured 
fractional logarithmic transmissions. 
20 Figure 31 shows the pulse height distribution of the photon spectrum 

transmitted through 15 cm muscle, 7 cm adipose and 1 cm bone simulating phantom. 
The spectrum was generated at 82 kVp, 0.15 mA and 750 second duration. 

Figure 32 (Phantom 3) shows the comparison of calculated values of the sum of 
products of the mass attenuation coefficients and areal densities and measured 
25 fractional logarithmic transmissions. 

Figure 33 shows the pulse height distribution of the photon spectrum 
transmitted through 15 cm muscle, 3 cm adipose and 2 cm bone simulating phantom. 
The spectrum was generated at 82 kVp, 0.144 mA and 750 second duration. 
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Figure 34 (Phantom 4) shows the comparison of calculated values of the sum of 
products of the mass attenuation coefficients and areal densities and measured 
fractional logarithmic transmissions. 

Figure 35 shows the pulse height distribution of the photon transmitted 
spectrum through 0.1 mm of copper (Cu). 

Figure 36 shows the pulse height distribution of the photons spectrum 
transmitted through 10 cm muscle, 3 cm adiopose and 1 cm bone simulating phantom. 
The spectrum was generated at 82 kVp, 0.029 mA and 1000 second duration. 

Figure 37 (Phantom 5) shows a comparison of calculated values of the sum of 
the mass attenuation coefficients and areal densities and measured fractional 
logarithmic transmissions. 

Figure 38 is a table showing the sum of areal densities of the phantoms 
calculated using Cramer's method. 

Figure 39 is a table showing calculated values of muscle areal densities. 

Figure 40 is a table showing calculated values of adipose tissue areal densities. 

Figure 41 is a table showing calculated values of bone areal densities. 

Figure 42 is a flow diagram illustrating a method for determining areal densities 
of bone, fat and soft tissue in a patient. 

DETAILED DESCRIPTION 

The present invention relates to a system and method for determining areal 
densities of bone, fat and muscle in a localized area of a patient. As illustrated in 
Figure 1, an apparatus can be used that includes a couch on which a patient lies, a 
scanning x-ray source and detector system and a computing system (not shown for 
purposes of clarity). Stored in or otherwise accessible to the computing system is a 
database or other collection of spectral curves that have been predetermined for the 
apparatus by experimental measurement, as described below. 

The coordinate system used in the following idea is rectangular and has +Y at 
the head of the patient lying supinely on a couch and -Y at the foot of the patient. +X 
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corresponds to the left side of the patient and -X corresponds to the right side of the 
patient. +Z is in the direction of the x-ray source above the patient and -Z is in the 
direction of the x-ray beam. 

The following method for determining areal densities of bone, fat and muscle in 
5 a localized area of a patient follows the flow diagram of Figure 42. After the database 
of curves has been created and stored as described below, the apparatus scans a 
localized area of a patient in a raster fashion, with the computer collecting data 
representing the measured intensity of the collimated x-ray beam after having passed 
through the patient at many points. At each point, the intensity of the beam is measured 

10 at a number of different filter thickness values interposed in the beam. An x-ray energy 
transmission curve is then computed from the collected intensity data. 

The x-ray source and detector system includes an x-ray machine fixedly 
coupled to an x-ray detector by a gantry (not shown in Figure 1 for purposes of clarity), 
and a scanning mechanism (also not shown) that moves the gantry laterally across the 

15 supine patient in either the ± x direction and then increments in the ± Y direction in a 
raster fashion. In this fashion, a collimated beam of x-rays is moved across a defined 
region of the patient that contains a bony structure. 

Suitable x-ray machines include those that can run continuously for several 
minutes in a kilo voltage range from approximately 40-50 kVp for infants and small 

20 children to 120-140 kVp for large and obese adults. The x-ray machine should be able 
to deliver up to 3,000 watts of electrical power, e.g., 20 mA at 140 kVp, for several 
minutes, e.g., 5-20 minutes. The x-ray beam is collimated to a small aperture which is 
aligned to hit a detector on the distal or back side of the patient. 

The detector system includes a rotating filter wheel or other suitable means for 

25 varying the intensity of the x-ray beam, suitable x-ray sensors and a collimator that 
directs the collimated beam upon the sensors. The filter wheel can be a metal (e.g., 
aluminum, copper, etc.) disk having sectors of different thicknesses. The thickness can 
increase from sector to sector such that the disk resembles a spiral staircase. 
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As described in further detail below, the x-ray transmission curve computed 
from the measured intensity data can be used to obtain a corresponding spectral curve, 
examples of which are shown in Figure 4. In the spectral curves, the number of x-ray 
photons (or spectrophotometer counts) at each energy present in the spectrum is plotted 
5 against the energy in keV. 

As noted above, an extensive database of pre-computed or predetermined 
spectral curves, which may resemble the curves shown in Figure 4, is accessible to the 
computer. These spectral curves or spectra have been pre-computed using a laboratory 
spectrometer in the system without a patient in place on the couch. In other words, the 

10 curves in the database describe the spectra of in-air beams, i.e., beams that pass only 
through air and are thus unattenuated. (The intensity of an in-air beam is referred to as 
Io in the equations described below.) The upper curve in Figure 2 illustrates the x-ray 
energy spectra for an in-air x-ray beam, and the lower curve illustrates the same 
spectrum after it has passed through the patient and has been attenuated. 

15 The computed x-ray spectra are stored in the database as a function of kVp and 

added filter. The added filters used in pre-computing the spectra need not be the 
spinning disk filter described above that is used when the system is actually in use with 
a patient. Rather, the added filters used in this pre-computing step can be variable 
filters of made of a suitable material. The variable filters can be composed of metal of 

20 different thicknesses or can be composed of different thicknesses of bone mineral, fat, 
and soft tissue. 

Figure 4 shows 11 x-ray spectra from 40 kVp to 140 kVp with 0.5 mm 
aluminum (Al) added filter. Figure 4 is intended merely to illustrate a number of 
spectra of the type that may be stored in the database; the database can actually include 
25 any of several thousand x-ray spectra. 

As noted above, a transmission curve computed from the measured intensity 
data can be used to obtain a corresponding spectral curve. Figure 5 illustrates this 
correspondence by showing transmission (attenuation) curves of each of the spectra 

197311 

8 



ATTORNEY DOCKET NO.: 21105.0005U2 

shown in Fig. 4 using an added aluminum (Al) filter. Note that this curve has a 
distinctive shape for each kVp. 

The upper range of added aluminum filter thickness in this example is 1 5 mm or 
1 .5 cm, but Al filter thicknesses of up several centimeters may be suitable. The typical 
5 adult body will have an attenuation corresponding to around 20 mm Al or two cm Al. 
The curve- matching process can thus start at around 15 to 20 mm Al. 

The database of spectral curves is searched until one is found that substantially 
matches the transmission curve defined by the collected data. As known in the art, an 
x-ray spectral curve relates beam energy to beam intensity at each beam energy. The 

10 bone density determination process can be started for any particular patient with a 
knowledge of the Io (e.g., the upper x-ray spectrum in Figure 2) as it is stored in the 
database and is determined from the kVp and added filter applied to the x-ray unit. 

As described above, when the system is in use on a patient, measurements are 
obtained with a spinning wheel Al filter of different thickness below the patient. 

15 Suppose the x-ray source is running at 140 kVp. A measured curve through 5.5 mm Al 
can be matched to the pre-computed curve for 140 kVp as a function of the filter by 
sliding it down the 140 kVp curve until it matches the original curve beginning at 6 mm 
Al and extending to 11.5 mm Al. This indicates that the patient or phantom at that 
point had an equivalent thickness of 6 mm Al. From this value, the computer looks up 

20 in the database the x-ray spectra at 140 kVp corresponding to 6 mm Al thickness and 
uses that for the lower x-ray spectrum in Figure 2. The computer then derives or reads 
from the database the x-ray spectrum from the group of x-ray spectra as shown in 
Figure 4 that is unique for that patient thickness and composition of soft tissue, fat and 
bone, at that kVp and added filter which enables the operations as described below. 

25 From the spectral curve found in the database, a curve is computed that relates 

the logarithm of the ratio of beam intensity in air divided by the measured beam 
intensity after passing through the patient to each of a plurality of photon energies in 
the beam. In other words, the computer is able to derive a spectrum such as the lower 
curve shown in Figure 2 and then calculate a curve such as that shown in Figure 3 
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below by knowing the original intensity at any photon energy Io(E) from the 
measurements and calculations we have already performed in the laboratory for that 
particular unit. Then, the computer computes the intensity 1(E) of the transmitted x-ray 
intensity through the patient for each photon energy. From the spectra, both 
transmission and unattenuated, the curve can be derived from log [Io(E)/I(E)j as shown 
in Figure #3. 

In the next step of the method, which is described in further detail below, the 
computer solves a matrix equation for the sum totals of the areal densities of bone, fat 
and muscle in response to predetermined mass attenuation coefficients for bone, fat and 
muscle at a group of three selected energies and in response to the logarithmic ratio 
obtained from the curve at each of the three selected energies. There are three energies 
in the group because that is the number needed for the matrix equation to determine 
three areal densities (bone, fat and muscle). The matrix equation is solved many times, 
each at a different group of three selected energies. 

Before the details of this step are described, it should be noted that the basis for 
determining the muscle, fat and bone mineral components at any x-y point or pixel is 
the following equation, which is well-known to persons skilled in the art to which the 
invention relates: 



Io(X,Y) is the x-ray intensity in air at the point or without a patient in place, I(X,Y) is 
the x-ray intensity in air with the patient in place, fi s ,f,b are the mass attenuation 
coefficients for muscle, fat and bone respectively, and X s , f ,b are the products of 
physical densities and thickness of muscle, fat and bone, respectively transverse by the 
photons beam. 

First, the following quantity is computed at each energy Ei in the x-ray spectra: 



I(X,Y) = I 0 (X,Y)exp(- Ms X s -ti f X f -fi b X b ) 



(1) 




= MbiX b +^ si X s +Li fi X f 



(2) 
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For the exemplary x-ray spectra shown in Figure 2, this energy range is 36.5 to 57.5 
keV. An exemplary curve described by this quantity at a number of energies Ei is 
shown plotted in Figure 3. The curve is 'the starting point for the mathematical solution 
5 to obtain the three component values. 

The exemplary curve shown in Figure 3 was obtained from 22 computations of 
Equation 2, but any suitable number of computations can be used. Note that in this 
example as many as 1,540 unique groups of three energy values along this curve can be 
selected to solve for the areal density of muscle, fat and bone mineral. (The number of 
10 permutations of three values selected from 22 values is 1,540.) It is believed suitable to 
use only 22 equations, with the values of photon energies separated by some optimal 
energy spacing. 

Three separate energies are arbitrarily selected from n energy intervals of the x- 
ray spectrum. Assuming narrow beam geometry, and that the object has three 
15 components, the three attenuation equations can be written for these energy intervals 
as: 

/, = I 0l exp{-ju b] X b -ju s] X 5 -ju fl X f ) 

I 3 =/ 03 exp(-ju b3 X b -Msi x s -M f 3 X /) 

where: Ioo,2,3) is the unattenuated intensity of the energy interval 1, 2 and 3, I(i,2,3) is the 
20 transmitted intensity in the energy interval 1, 2 and 3. The coefficients, |i s ,f;b , are the 
mass attenuation coefficients in cm /g of each component whose surface densities are X 
in g/cm 2 , and b, s, f are indices for bone, soft tissue and fat respectively. The equations 
can be written in matrix notation as: 



Mm 




/"/, 




'x b ' 




"ln(/ 01 //,)- 




M s 2 


Mf2 


X 






ln(/ 02 // 2 ) 


y"63 


Vs3 


M f3 _ 




x f 




ln(/ 03 // 3 ) 



25 or in the simpler form A x X=B (5) 
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and the solution can be obtained using the inverse 
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5 

The linear equation contains three unknowns: X b , X s , and Xf. Several mathematical 
methods have been applied to solve such a linear equation. 

The determinant of matrix A in Equations 4 and 5 above has a very small value. 
This makes it very difficult to solve as it is ill conditioned. If perfect input data were 

10 put into matrix B then a perfect solution would be obtained, but this is not possible for 
any kind of measured data. The individual values of the components obtained by this 
method will vary widely for each solution and they will range from large negative 
values to large positive values for the different components. These values make no 
physical sense in many cases. By studying the solutions, it was discovered in 

15 accordance with the present invention that even though the individual solutions were 
not correct, the sum of the individual solutions for the three components was almost 
correct. ("Almost correct" in the sense that it varied at most only a few percent or 
always less than 5% from the known correct sum solution for phantoms that were 
measured in the laboratory.) The novel method of the present invention for 

20 determining the three components is based upon this discovery. For example, two of 
the values can be negative and one can be positive, but it was discovered that their sum 
would be very close to the correct and true sum. This was determined from phantom 
measurements with accurately known chemical compositions. 

The computer, operating in accordance with suitable programming to effect the 

25 method described herein, takes the sum of the individual solutions and estimates the 
areal density of the individual components. The starting point for this process is patient 
parameter-dependent. For example, a first estimate can be made for the fraction of 
muscle, fat and bone mineral along a ray which is dependent upon the patient's sex, 
age, level of obesity and exogenous factors such as known osteoporosis, diseases, etc. 
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An exemplary software code excerpt illustrating the iterative process is shown in Fig. 6. 
In the code, Kgl represents the estimated fraction of total areal density for muscle, Kg2 
represents the estimated fraction of areal density for fat or adipose tissue, kg3 
represents the estimated fraction for bone, sumrhoxg denotes the average sum of the 
5 different determinant or matrix solutions or correct value, p5 denotes a particular 
phantom or patient, and 1,2,3 in the array denote soft tissue, fat and bone, respectively. 
Sumdifbest is initialized to a large value to start the loop. The loop iterates to 
determine the best fraction of total areal density for soft tissue, fat and bone. 
Conceptually, on each iteration, the computer is generating a curve in response to 
10 estimates of areal densities of bone, fat and muscle and then comparing the generated 
curve to the logarithmic curve until a match within predetermined threshold limits is 
identified. Persons skilled in the art to which the invention relates will readily be 
capable of providing a computer with suitable programming to effect the above- 
described steps. 

15 As in dual x-ray bone densitometry, the x-ray beam has passed through a 

thickness of tissue adjacent to a bone structure where the x-ray beam does not actually 
encounter any bone structures, i.e. right adjacent to the spine in the body below the 
diaphragm and above the pelvis. In such a case, there are only two components: soft 
tissue and fat. If the transmitted x-ray spectrum is known at this point, the low and 

20 high energy can be selected from the x-ray spectrum and solved to obtain the area 
density of soft tissue as in dual x-ray absorptiometry as: 

x = InQoL 1 h ) ~ [hifloH 1 r H ) W I l^fH (7) 

25 The subscripts L and H correspond to lower and higher photon energies. 
Similar equation can be written for areal density of adipose tissue. 

As the scanning x-ray beam passes into the region containing bone, the solution 
values for muscle and fat outside the bone can be taken and subtracted from the sum 
values with bone included to obtain a first estimate for the bone mineral values. The 
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values of areal density obtained for each component are passed into a an iteration 
routine whereby the type of curve shown in Figure 3 is calculated and compared to a 
measured such as the measured curve shown in Figure 3 curve. The sum of the 
absolute difference at point by point energy comparison is calculated. Then, in an 
5 iterative fashion, the measured fractional areal densities of the different components 
can be varied slightly in both the plus and minus directions and a new curve is 
calculated. If the sum of the absolute difference of this curve is greater than the 
previous absolute difference, the areal densities can be perturbed in the opposite 
direction. In this fashion, each fractional areal density can be perturbed to finally 
10 obtain those fractional values of the areal density that yield the least sum of the 
absolute difference, energy point by energy point, of the calculated and measured 
curves. 

In cases in which values for muscle and fat cannot be obtained, i.e. the beam 
does not go through a bone- free region adjacent to bone, a number of matrix solutions 
15 (e.g., from 1 to several thousand) can be computed and an average value of their sum 
computed. 

While the description set forth above is believed sufficient to enable persons 
skilled in the art to which the invention relates to make and use the invention, the 
following discussion of the theoretical and practical underpinnings of the invention 

20 may be of interest to the reader. 

A specific underpinning of the invention involved a study with experiments in 
which the x-ray spectra incident and transmitted through a phantom simulating tissue 
were measured using a semiconductor detector and multichannel analyzer (not shown). 
The measured spectra were used to calculate the mass attenuation coefficients of the 

25 phantoms simulating tissue as a function of energy and compare them with theoretically 
derived values. 

Another specific underpinning of the invention involved experiments in which the 
x-ray spectra incident and transmitted through the multiple tissue simulating 
compositions was measured. Spectra were obtained for a variety of phantom 
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compositions. The fractional logarithmic transmissions were calculated for each 
phantom and compared with the sum of the products of mass attenuation coefficients 
and areal densities of the phantom components. 

Still another specific underpinning of the invention involved developing a 
5 computer program to calculate the areal densities of the phantom components (muscle, 
bone, adipose tissue) from the measured spectra. 

The spectroscopy system used in these studies (not shown) employed a planar, 
low energy, high purity germanium (LEGe) detector cooled in liquid nitrogen, 
Canberra model GL 0510. The detector had an active area of 500 mm 2 , diameter of 

10 25.2 mm and a thickness of 10 mm. The detector was equipped with a 0.15 mm 
beryllium (Be) window and the distance from the beryllium window to the germanium 
crystal was 5 mm. The bias voltage of the detector was 1500 volts. The preamplifier 
was built into the detector assembly. The spectroscopy system used a Canberra 
amplifier model 2022 and a Canberra high voltage supply model 1000. The output 

15 from the amplifier was connected to the multichannel analyzer, such as an Ortec- 
Norland 5600 model having a 100 MHz ADC. 

The energy calibration of the multichannel analyzer was performed by exposing 
the detector to known energies of Co-57 (14.41 keV), Am-241 (59.536 keV) and Ba- 
133 (80.99 keV) sources. The photopeaks were centered in the energy channel using 

20 the built in multichannel analyzer function. The conversion gain of the multichannel 
analyzer was set to 1024 and the amplifier was set to display 1 keV per channel, i.e. 
channel 20 measured photons from 19 to 20 keV. The linearity of the multichannel 
analyzer (the relationship between gamma-ray energy and multichannel analyzer 
channel number) was measured with the same sources as used for the energy 

25 calibration. 

The data acquired in this study was measured using a constant potential dental 
x-ray unit (not shown), Keystone model 2497-1 00 Y, equipped with a cone. The system 
operated in the 60-90 kVp range with tube currents of 0.01 mA to 1 mA. There is a 
very small ripple associated with the constant potential output (< 3%). The unit 
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employs a thick tungsten target with an anode target angle of 15 degrees. The beam is 
modified by 3 mm aluminum total filtration (1 mm inherent and 2 mm added). The 
kVp and tube current are displayed during the exposure at the control panel. During the 
data collection, the x-ray target-detector distance was 130 cm. To reduce the primary 
photon fluence rate, the x-ray beam was collimated in three different steps. The source 
was equipped with a lead collimator with a diameter of 15 mm and thickness of 40 mm. 
A lead collimator with an aperture diameter of 1 mm and lead thickness of 5 mm plus 2 
mm of aluminum was used at the detector so that photon penetration through the 
collimator wall was negligible for operating tube potentials. Additionally, a lead 
collimator with the diameter 2 mm and thickness of 65 mm was positioned 66 cm from 
the target. 

An important element in this study was the stable and reproducible energy 
spectrum of the x-ray unit. To achieve a sufficient number of photons without 
overheating the x-ray unit, the exposure time was set at 50 seconds. The multichannel 
analyzer operated in the "add all" mode and added photons emitted in the successive 
exposures. To check for the reproducibility of x-ray spectrum, ten exposures of 50 
seconds duration were taken at 85 kVp and 0.01 mA current. These are the x-ray 
spectra shown in Figure 13. To reduce the dead time, the spectrum was measured with 
7.0 cm of acrylic in the beam. The dead time as displayed by the multichannel analyzer 
was 0-1%. The reproducibility was obtained by calculating coefficients of variation for 
the selected energy channels. The standard deviation of the ten 50 second exposure 
spectra was calculated from the formula: 



Where: SDy is the standard deviation of the mean value; Ij is the observed 
number of counts in the selected energy interval and N is the number of exposures. 
The coefficient of deviation for the selected energy intervals is: 



10 



SD 7 = 



(8) 



N 




(9) 
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As mentioned above, for the measurements of the mass attenuation coefficient 
as well as the bone mineral measurements, the photons collected for several 50 second 
exposures were added. The standard deviation of the mean and coefficient of variation 
for six 500 seconds exposures (which is the sum often 50 seconds exposures) at each 
photon energy were also calculated. 

Tissue simulating phantoms were manufactured by Computerized Imaging 
Reference Systems, Inc. Figure 7 shows the composition and fraction by weight of the 
constituent elements provided by the manufacturer for the three types of tissue used in 
this study. The mass attenuation coefficients of the tissue simulating material 
components was taken from Hubbell and Seltzer and interpolated at 1 keV intervals 
using the Newton interpolation divided difference formula. {See Hubbell, et ah "Tables 
of X-ray mass attenuation coefficients and mass energy absorption coefficients 1 keV 
to 20 MeV for elements z=l to 92 and 48 additional substances of dosimetric interest." 
NIST, 1995.) The table values were derived at mid-channel values, i.e. 19.5 keV for 
channel 20. The values of the mass attenuation coefficients of the tissue simulating 
components were calculated by adding the values of the mass attenuation coefficients 
of the elements according to their proportions by weight. 

E iE) = ^J?-{E) 
P i VP Ji 



(10) 



where: w\ is the fraction by weight of the i-th element and (|^/p)i are the interpolated 
values of atomic element of the mass attenuation coefficients at the photon energy E. 

The table of Figure 7 shows the densities and fractional composition by weight 
of tissues included in the measurements. 

The response function of the spectrometer is the analyzer pulse height 
distribution. To obtain an accurate determination of the incident photon spectrum, the 
measured pulse height distribution needs to be corrected for detector distortion. The 
measured values used in the experiment are the pulse height distribution not corrected 
for detector characteristics. This was based on assumption, that both primary and 
transmitted spectra were distorted in the same way by the measurements. The 

197311 

17 



ATTORNEY DOCKET NO.: 21105.0005U2 

distortion will cancel out in both denominator and numerator. We believe that this 
assumption is correct if both measurements were made with the same dead time. 

The materials simulating the muscle, adipose tissue and bone were placed 
behind the middle collimator, 66 cm away from the tube target. To achieve the same 
dead time as displayed on multichannel analyzer all measurements were taken with 
7 cm of acrylic in the primary beam. The acrylic was placed behind the middle 
collimator. The transmitted beam was measured by placing the tissue simulating 
materials behind the acrylic. Exposures of 500 seconds duration were taken for both 
the primary and transmitted spectra. For all measurements, the dead time, as displayed 
on the multichannel analyzer, was 0-1%. The mass attenuation coefficients were 
calculated knowing the overall thickness and physical density of the tissue from the 
attenuation equation: 



where: Io is the number of photons measured in the energy channel with only the 
acrylic in the beam and I is the number of photons measured with acrylic and muscle, 
adipose, or bone simulating phantoms, p is the physical density, in g/cm 3 , and x is the 
thickness of the tissue transversed by the photons, in cm. The measured mass 
attenuation coefficients were fitted to a second order polynomial. 

Assuming a perfect detector, negligible error in the thickness measurement and 
ignoring any effect due to scatter radiation, the uncertainty associated with the 
measurement of the mass attenuation coefficients was calculated using propagation of 
errors from the formula: 



where: a^ 0 ) is the standard deviation, p is the physical density of the material being 
measured,Iois number of counts per channel for the primary beam and I is the number 
of counts for the transmitted beam. 



//_ln(/ 0 /7) 



(11) 



P P* 




(12) 
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The shape of the spectrum, when the x-ray beam was attenuated by 15 cm of 
water is shown in Figure 9 and the primary spectrum is shown in Figure 8. On 
comparison with the primary spectrum the maximum is shifted to approximately 
57 keV and photons below 30 keV are almost completely absorbed. This is the reason 
5 that the data analysis was restricted to the energy range of 36.5-57.5 keV on the x-ray 
spectrum. 

Figures 8 and 9 illustrate spectra used to determine the energy range of interest. 
Figure 8 shows the pulse height distribution of the photons of the primary spectrum for 
87 kVp, 0.01 mA tube current and 500 second duration obtained with the acrylic water 
10 container and an additional 3.2 cm of acrylic in the x-ray beam. Figure 9 shows the 
pulse height distribution of the photon spectrum obtained after filling the container with 
1 5 cm of water. 

Five different phantoms were assembled from the three components simulating 
human tissue. The geometry consisted of stacks of three slabs with varying thickness 

15 of each of three materials with the adjacent slabs in contact with each other. In four of 
the phantoms, bone mineral areal densities were 1.084 g/cm 2 . The value of about 
1 /cm 2 corresponds to the value typically encountered in the measurement of the spine 
and femoral neck. {See Dunn, et aL 9 "Measurement of bone mineral content in human 
vertebrae and hip by dual photon absortiometry." Radiology 156, 203 (1985); Wahner, 

20 et al 9 "Dual-photon Gd-153 absortiometry of bone." Radiology 156, 203 (1985). 
Measurements of bone mineral density over a range of soft tissue compositions were 
made by changing the relative thickness of muscle and fat. The simulated muscle had 
thicknesses of 10, 15 and 20 cm. The adipose tissue mass percentage of soft tissue in 
the assembled phantoms ranged from 13% to 40%. These values span the normal range 

25 of fat measured in the body. {See Wang, et aL 9 "Body fat from body density: 
underwater weighting vs. dual-photon absortiometry." Am. J. Physiol. 256, 1989). The 
composition and thickness of phantom components is shown in Table 2. To prevent 
saturation of the spectrometer system and keep the detector pile-up effect at a 
negligible level the open beam (Io) was measured with the lowest tube current available 
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10 



15 



20 



25 



(0.01 mA) and with 0.1mm of Cu in the beam. The data obtained was corrected for the 
copper attenuation and the ratio of currents. Phantom measurements were obtained 
with higher tube currents. During the phantom measurements, the tube current varied 
from 0.029 mA to 0.343 mA. 

By knowing the mass attenuation coefficients and areal densities of the three 
phantom components, the sum of the products of mass attenuation coefficients and 
areal densities of the phantom components can be calculated from formula 13 (right 
side of equation). 



In 



x 0i 



= ^bi X b+Msi X s+^fiX f 



(13) 



where I 0 i corresponds to the number of photons in the primary beam at energy interval i, 
Ii is the transmitted number of photons through the phantom at the same energy 
interval, /x b , s ,f is the mass attenuation coefficients for bone, muscle and adipose tissue, 
respectively, and X b , s ,f are the areal densities of the phantoms components. The table of 
Figure 10 shows the values of areal densities and the compositions of the phantoms. 

From the measured spectra, the fractional logarithmic transmission (left side of 
the equation 13) was calculated for the energy range of 36.5-57.5 keV and compared 
with the sum of the products of mass attenuation coefficients and areal density of 
phantom components (right side of the equation 13). 

Three separate energies can be randomly selected from n energy intervals of the 
x-ray spectrum. Assuming narrow beam geometry, and that the object has three 
components, three attenuation equations can be written for these energy intervals as: 
/, =/ 01 exp(-/i bl X b -fi sl X s -f4 fl X f ) 

I 2 =I 02 exp(-M b2 Xt> -M s2 X s -Mf 2 X f ) 04) 

where: Io(i,2,3) is the unattenuated intensity of the energy interval 1, 2 and 3, I(i,2,3) is the 
transmitted intensity in the energy interval 1, 2 and 3. The coefficients, n St f f b > are the 
mass attenuation coefficients in cm 2 /g of each component whose surface densities are X 
in g/cm 2 and b, s, f are indices for bone, muscle and fat respectively. 
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The equations can be written in matrix notation as: 
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or in the simpler form A x X=B 

(16) 

5 and the solution can be obtained using the inverse [/z]~\ 
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The linear equation contains three unknowns: X b , X s , and Xf. Several mathematical 
methods have been applied to solve this linear equation. 

As previously stated, the energy range selected on the x-ray spectra for use in 
10 this project contains twenty two channels in the energy range of 36.5 to 57.5 keV. 

To solve the linear equation it is necessary to establish that the determinant of 
the matrix A, which stores the mass attenuation coefficients, is different from zero. 
Because the difference in the mass attenuation coefficients is small the matrix is ill 
conditioned and the determinant of the matrix as reported by several scientists is very 
15 small. (See Lehman, et al. 9 "Generalized image combinations in dual KVp digital 
radiography." Med. Phys. 8. 1981 (659-667); Kotzki, et al, "Theoretical and 
experimental limits of triple photon energy absorptiometry in the measurement of bone 
mineral." Phys. Med. Biol 36(4), 429-437, 1991; Sutcliffe, "A review of in vivo 
experimental methods to determine the composition of the human body." Phys. Med. 
20 Biol. 41, 1996 (791-833).) 

For example, when calculating the determinant of the matrix at 1 keV interval, 
matrix A, which stores the mass attenuation coefficients for bone, muscle and adipose 
tissue, respectively, at the energies 36.5, 37.5 and 38.5 keV is: 
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0.3799 0.2891 0.2505 
0.3668 0.2821 0.2467 
0.3543 0.2754 0.2429 

The determinant of the matrix is -1. 58427- 10~ 7 . 

Increasing the energy interval to 6 keV, which corresponds to mass attenuation 
coefficients values at the energies 36.5, 42.5 and 48.5 keV matrix A is now: 

0.3799 0.2891 0.2505" 
0.3104 0.2516 0.2293 
0.2631 0.2250 0.2129^ 

The calculated determinant of the matrix shown above is -1.437997 10" 6 . 

10 A program was written in True Basic using an IBM computer to calculate the 

areal densities of the phantom components, matrix X. Nevertheless, persons skilled in 
the art to which the invention relates will readily be capable of writing or otherwise 
providing such a program in other languages or for use on other computer systems. 
The program employs a guessing routine to minimize the difference between guess 

15 values and measured values of the fractional logarithmic transmissions. The fractional 
logarithmic transmissions through the phantoms of three components that were 
measured as described above in accordance with another specific underpinning of the 
invention, were used as the input [B] matrix. Matrix [A] is the two dimensional matrix 
representing the mass attenuation coefficients of phantom components. The program 

20 allows the user to select the energy interval. Starting with the lowest energy in the 
energy range of interest (36.5 keV), the next energy is 36.5 keV plus the selected 
energy interval, and the third is 36.5 keV plus two intervals. When the energy reaches 
the highest value of the selected energy range (57.5 keV) the loop routine takes it back 
to the lowest value of energy until the twenty two solutions are calculated. For 

25 example, when the selected energy interval is 4 keV, the first row in matrix A stores 
mass attenuation coefficients for bone, muscle and adipose tissue at the energies of 36.5 
keV, second row at the energies of 40.5 keV and third at the energy of 44.5 keV. In the 
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second set of equations, the first row stores the values of mass attenuation coefficients 
at the energy of 40.5 keV, the second row at the energy of 44.5 and third at the energy 
of 48.5 keV. In this program, twenty two sets of equations were calculated. 

The sets of twenty two equations are first solved using the Cramer's rule and the 
5 calculated twenty-two values of areal densities of the phantom components are 
averaged. When solving the equations using Cramer's method, some solutions are 
obtained which do not make physical sense, as some of the calculated areal densities 
have negative values. Although, the calculated average values of areal density of the 
phantom components are sometimes negative their sum is always very close to the 

10 calculated values of the sum of the areal densities of the phantoms. The real value of 
the sum of the phantom areal densities was calculated as Xb + X s + Xf. Phantom 1 
contains 20 cm of muscle simulating component, 1 cm of bone and 3 cm of adipose 
tissue. This corresponds to 21.24 g/cm 2 of areal density of muscle; 1.084 g/cm 2 of bone 
and 2.901 g/cm 2 of adipose tissue. When solving a set of linear equations using 

15 Cramer's rule and an energy interval of 4 keV, the calculated values of areal densities 
for the first set of equations are +20.411 g/cm 2 , -0.542 g/cm 2 , and +5.837 g/cm 2 for 
muscle, bone and adipose tissue, respectively. The sum of the calculated areal density 
using the Cramer's method is 25.706 and the calculated true value is 25.225 g/cm 2 . 
Phantom 5 contains 10.62 g/cm 2 of areal density of muscle, 1.084 and 2.901g/cm 2 of 

20 bone and adipose tissue, respectively. The calculated values of areal density using the 
same, 4 keV energy intervals for the first set of equations and the Cramer's method are 
13.243, 0.224 and 1.176 g/cm 2 for muscle, bone and fat, respectively. The sum of the 
calculated areal density is 14.643 g/cm 2 with true value being 14.605 g/cm 2 . 

For every phantom, the fraction of areal density of the phantom components as 

25 a fraction of the sum of areal densities of muscle, adipose tissue and bone was 
calculated as: 

f s = — (18) 

X 5 +X b +X f 
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and the fraction of adipose tissue can be calculated as ff =l-(f s + fb). 
where: f s , fb and ff are the fractions of muscle, bone and adipose tissue, respectively. X s , 
Xb and Xf are the areal densities of the muscle, bone and adipose tissue, respectively. 
5 Areal density, as mentioned previously, is the product of the physical density and 
thickness of the phantom component transverse by the photons in g/cm 2 . 

For the five phantoms used in this project, the muscle fraction ranges from 
0.669 to 0.842, the adipose tissue fraction ranges from 0.114 to 0.286 and the bone 
fraction ranges from 0.043 to 0.103. For example, the areal density of muscle in 

10 phantom 1 is 21.24 g/cm 2 , and the sum of areal density of all phantom components is 
25.225 g/cm 2 . This corresponds to the fraction of muscle in this phantom of 0.842. 
The same phantom contains 1.084 g/cm 2 of bone, which corresponds to the fraction of 
bone in this phantom of 0.043. The fraction of the adipose tissue is calculated by 
subtracting from one the sum of fractions of muscle and bone fractions. 

1 5 For each phantom, the program assigns the range of the fraction of areal density 

of the phantom components by adding and subtracting 0.01 from the calculated 
fraction. For example the fraction of the bone in phantom 1 is 0.043. The calculated 
range varies form 0.033 to 0.053. Knowing that fraction of 0.043 corresponds to 1.084 
g/cm 2 , the fraction of 0.033 corresponds to the areal density of the bone of 0.832 g/cm 2 

20 and fraction of 0.053 to 1.336 g/cm 2 . The range of 0.7-1.4 g/cm 2 of bone mineral 
density covers the normal spine and femoral neck scanning. The change of the fraction 
by only 0.01 in both directions change the fraction of the bone by 23% in both minus 
and plus direction. 

In the calculation loop of the computer program, the first estimate is calculated 
25 as the sum of the products of the fractions of the phantom components and calculated 
sum of the areal densities and mass attenuation coefficient of the phantom components 
at selected energies. The value of the fraction is perturbed in 0.001 step in both minus 
and plus directions within calculated fraction range to achieve the minimum difference 
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between the calculated sum of the products of the areal densities and mass attenuation 
coefficients and measured fractional logarithmic transmission. Measured logarithmic 
transmission for each phantom at the selected energy range (36.5-57.5 keV) is matched 
to the calculated sum of the products of the fractions of phantom components, mass 
5 attenuation coefficients of the components and sum of areal density for that phantom. 
The program assigns the best fraction for each of the phantom component for which the 
difference between the sum of the products of these fractions and mass attenuation 
coefficients of this components and sum of the areal density of the phantom and 
measured fractional logarithmic transmission is the smallest. 

10 The areal density of the phantom components is calculated as the product of the 

best fraction of that component and the sum of areal densities. 

The results are described below. The calibration spectra of Co-57 and Am-241, 
plotted on the same graph, are shown in Figure 11. The channel numbers for Co-57, 
Am-241 and Ba-133 sources plotted against photon energy in keV are shown in 

15 Figure 12. 

Ten of the x-ray spectra obtained in 50 second exposures are shown in 
Figure 13. The coefficients of variation were calculated for selected energies in order 
to determine the variability of the mean. The mean number of counts, standard 
deviation and coefficient of variation are shown in the table of Figure 14. The average 
20 value of coefficient of variation calculated from the formula 1 is 0.0109. 

Ten of the 500 seconds spectra, which are ten 50 seconds exposures added are 
shown in Fig. 15. The table of Figure 16 shows the standard deviation, coefficients of 
variation of the 500 seconds exposures. The average value of the calculated coefficient 
of variation is 0.0039. 

25 Figure 17 shows the theoretically derived values of the mass attenuation 

coefficients for muscle, adipose tissue and bone simulating phantoms used in this study 
plotted as the function of energy in the energy range of interest. 

The values of the components of mass attenuation coefficients were taken from 
Hubbell and interpolated at the middle of the channel at 1 keV intervals. The values of 
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mass attenuation coefficients of tissues calculated from the formula 10 and plotted as 
the function of energy in the energy range of interest are shown in Figure 17. 

Figure 18 shows the primary and transmitted spectra through 5 cm of muscle 
simulating phantom. In Figure 19 the experimentally obtained mass attenuation 
coefficients are compared with theoretically determined values. The measured mass 
attenuation values were fitted to a second degree polynomial. 

The average absolute percent difference between the measured and the table 
values of the mass attenuation coefficients at 1 keV intervals were calculated from the 
formula below: 

r 



abs 

Average % difference = ^ 
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The average percent difference between the experimentally and theoretically 
determined values calculated from the formula 20 in the energy range of 36.5-57.5 is 
0.91%. The values of experimentally and theoretically derived values of mass 
attenuation coefficients in the energy range of interest, ratios of experimentally and 
theoretically derived values and uncertainty in the mass attenuation coefficient per 
energy channel are shown in the table of Figure 20. Using Formula 12, the average 
uncertainty in the mass attenuation coefficient of the muscle simulating phantom is 
0.0037 cm 2 /g. 

The primary and transmitted pulse height distributions through the 3 cm of the 
adipose tissue simulating phantom are shown in Figure 21 . The measured, theoretically 
derived and second degree polynomial fitted to the measured data is shown in 
Figure 22. The calculated average percent difference between the experimentally and 
theoretically determined values calculated from equation 12 was 0.54 %. The average 
calculated uncertainty in the experimentally determined values was 0.0056 cm 2 /g. The 
table of Figure 23 shows the values of experimentally and theoretically derived values 
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of mass attenuation coefficients in the energy range of interest, ratios of experimentally 
and theoretically derived values and uncertainty in the mass attenuation coefficient per 
energy channel. 

The primary and transmitted pulse height distributions through 2 cm of the bone 
5 simulating tissue are shown in Figure 24. Comparison of measured, theoretically 
derived, and second degree polynomial fitted to the measured data is shown in 
Figure 25. The calculated average percent difference between the experimentally and 
theoretically determined values was 0.81 %. The average calculated uncertainty in the 
experimentally determined values was 0.0069 cm 2 /g. The table of Figure 26 shows the 
10 values of experimentally and theoretically derived values of mass attenuation 
coefficients in the energy range of 36.5-57.5 keV, ratios of experimentally and 
theoretically derived values and uncertainty in the mass attenuation coefficient per 
energy channel. 

The following results were obtained for measurements through three component 
1 5 phantoms and calculation of the fractional logarithmic transmissions. 

Bone 100 (Computerized Imaging Reference Systems, Inc.) was used in all of 
the phantoms assembled from the tissues simulating materials. Bone 100 corresponds 
to 100 mg/cm 3 of hydroxyapatite (Cai 0 (PO 4 )6(OH) 2 ). A range of hydroxyapatite of 
trabecular bone from 50-200 mg/cm 3 is typical for spinal bone mineral measurements. 
20 Phantom 1 contains 20 cm of the muscle simulating material, 3 cm of adipose 

tissue and 1 cm of bone. This corresponds to areal densities of 21.24 g/cm 2 of muscle, 
2.901 g/cm 2 of adipose tissue and 1.084 g/cm 2 of bone, respectively. The overall 
thickness of the phantom was 24 cm. The spectrum transmitted through the phantom is 
shown in Figure 27. Comparison of the calculated values of sum of the products of 
25 areal densities and mass attenuation coefficients (right side of equation 13) and 
measured values of fractional logarithmic transmission (left side of equation 13) is 
shown in Figure 28. The average percent difference calculated from formula 20 at the 
energy range of 36.5-57.5 was 0.63%. 
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To examine the influence of body thickness on the measured bone mineral 
value, the second phantom contained 15 cm of muscle and the same amount of bone 
and adipose tissue as the previous phantom. This corresponds to areal densities of 
15.93 g/cm 2 of muscle, 2.901 g/cm 2 of adipose tissue and 1.084 g/cm 2 of bone. The 
5 spectrum transmitted through the phantom is shown in Figure 29. The comparison 
between the calculated sum of the products of the mass attenuation coefficients and 
areal densities (left side of equation 13) of the three components phantom and 
fractional logarithmic transmission (right side of equation 13) is shown in Figure 30. 
Calculated from formula 20, the average percentage difference between the calculated 

10 and measured value is 0.26%. 

Figure 30 shows the comparison of calculated values of the sum of products of 
the mass attenuation coefficients and areal densities for phantom 2 and measured 
fractional logarithmic transmissions. 

Phantom three contains 15 cm of muscle, 7 cm of adipose tissue and 1 cm of 

15 bone. This corresponds to areal densities of 15.93 g/cm 2 of muscle, 6.769 g/cm 2 of 
adipose tissue and 1 .084 g/cm 2 of bone. The pulse height distribution of the photons 
spectrum transmitted through the phantom is shown in Figure 31. Calculated from 
equation 13, the sum of the products of the mass attenuation coefficients and areal 
densities of the phantom components (left side of equation 13) are compared with the 

20 logarithmic transmission factor in Figure 32. The average percent difference calculated 
from Formula 13, over the energy range of interest is 0.47 %. 

Phantom four contains 15 cm of muscle, 3 cm of adipose tissue and 2 cm of 
bone, which corresponds to the areal densities of 15.93 g/cm 2 for muscle, 2.901 g/cm 2 
for adipose tissue and 2.168 g/cm 2 for bone. The pulse distribution of photons 

25 transmitted spectrum through the phantom is shown in Figure 33. The comparison of 
the calculated sum of the products of the mass attenuation coefficients and areal 
densities of phantom components and fractional logarithmic transmission is shown in 
Figure 34. Calculated from the equation 20, the average difference between measured 
and calculated values is 0.21%. 
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For the four phantoms described above the primary spectrum was obtained at 
the lowest tube current of 0.01 mA, at the same kVp as the phantoms measurements. 
To reduce photon fluence the pulse height photon spectrum was measured with 0.1 mm 
of copper (Cu) in the beam. The pulse height distribution of photons transmitted 
5 through 0.1 mm of Cu is shown in Figure 35. The measured pulse height photon 
spectrum was corrected for 0. 1 mm Cu attenuation and the ratio of tube current for each 
phantom. Tube currents for the phantom measurement were increased to 0.343 mA for 
phantom 1, 0.152 mA for phantom 2, 0.15 mA for phantom 3 and 0.144 mA for 
phantom 4. The dead time for the phantom measurements and the primary beam 

10 displayed on the multichannel analyzer was 0-1%. An exposure of 750 seconds was 
performed for each measurement. 

Figure 35 shows the pulse height distribution of photons transmitted spectrum 
through 0.1 mm of copper (Cu). 

Phantom five contains 10 cm of the muscle, 3 cm of adipose and 1 cm of bone. 

15 This corresponds to areal densities of 10.62 g/cm 2 of muscle, 2.901 g/cm 2 of adipose 
tissue and 1.084 g/cm 2 of bone. The pulse height distribution of the photon spectrum 
transmitted through the three components phantom is shown in Figure 36. The 
comparison of the calculated sum of the products of mass attenuation coefficients and 
the areal densities of phantom components and the fractional logarithmic transmission 

20 is shown in Figure 37. Calculated from equation 20, the average difference between 
calculated and measured values in energy range of 36.5-57.5 keV was 0.36%. The 
measurement of the primary beam for this phantom was performed with a tube current 
of 0.01 mA 82 kVp and 7 cm of acrylic in the beam. During the transmitted beam 
measurements 7 cm of acrylic was added to the beam. The exposure time 1000 second 

25 duration was used for both primary and phantom spectra. During the phantom 
measurement, the tube current was increased to 0.029 mA. 

Figure 36 shows the pulse height distribution of the photons spectrum 
transmitted through 1 0 cm soft tissue, 3 cm adipose and 1 cm bone simulating phantom. 
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The spectrum was generated at 82 kVp, 0.029 mA and 1000 second duration. 
Additionally, 7 cm of the acrylic in the x-ray beam. 

Figure 37 shows a comparison of calculated values of the sum of the mass 
attenuation coefficients and areal densities for phantom 5 and measured fractional 
5 logarithmic transmissions. 

When calculating the sum of areal density of the phantoms using the Cramer's 
method, as previously stated, the sum of the calculated areal density is close to the sum 
of true values. The table of Figure 38 shows the sum of the areal densities of the 
phantoms, and the calculated average values using the Cramer's method. The values 

10 were obtained by increasing the energy intervals from 1 keV to 10 keV. The Cramer's 
method overestimates the sum of the areal densities for all the phantoms. Increasing the 
energy interval from 1 keV to 10 keV the percent difference between the calculated 
average sum of the areal densities and the true values of sum of the areal densities is 
2.75% for phantom one, 0.21% for phantom two, 3.05% for phantom three, 1.7% for 

15 phantom four and 1.75% for phantom five. The standard deviation of the mean values 
is 0.082, 0.028, 0.154, 0.052 and 0.058 g/cm 2 for phantom 1, 2, 3, 4 and 5, respectively. 

Using the computer program and increasing the energy interval from 1 to 10 
keV the areal densities and the average values for muscle, adipose tissue and bone were 
calculated for each phantom and are shown in Figures 39, 40 and 41, respectively. 

20 Increasing the energy interval from 1 to 1 0 keV the average values of the areal densities 
of the muscle is overestimated by 2.26% for phantom 1, 0.1 1% for phantom 2, 3.18 % 
for phantom three, 1.57% for phantom four and 1.17% for phantom five. The 
calculated standard deviation of the average value is 0.069, 0.021, 0.089, 0.047 and 
0.030 g/cm 2 for phantoms 1, 2, 3, 4 and 5 respectively. 

25 The mean value of the calculated areal densities of the adipose tissue, shown in 

Figure 40, is overestimated by 7.21% for phantom 1, 3.20% for phantom 2, 5.11% for 
phantom 3, 4.98 % for phantom 4 and 6.12% for phantom 5. The standard deviation of 
the calculated mean is 0.010 g/cm 2 for phantom 1, 0.023 g/cm 2 for phantom two, 0.058, 
0.016 and 0.038 g/cm 2 for phantoms three, four and five respectively. 
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Increasing the energy intervals from 1 to 10 keV the calculated values of the 
areal densities of the bone mineral and the mean values for each phantom are shown in 
Figure 41. The calculated average value of the areal density of the bone is higher by 
0.42% a true value for the phantom 1. The average calculated value of the bone 
5 mineral density is underestimated by 6.31%, 2.71%, 1.71% and 4.21% for phantoms 2, 
3, 4 and 5, respectively. The calculated standard deviation of the mean value is 0.003, 
0.017, 0.017, 0.012 and 0.010 g/cm 2 for phantoms 1, 2, 3, 4 and 5, respectively. 

The table of Figure 38 shows the sum of areal densities of the phantoms 
calculated using Cramer's method. 
10 The table of Figure 39 shows calculated values of muscle areal densities. 

The table of Figure 40 shows the calculated values of adipose tissue areal 
densities. 

The table of Figure 41 shows the calculated values of bone areal densities. 

For the accurate determination of bone mineral density using the dual or triple 

1 5 energy method, the determination of mass attenuation coefficient of body components 
is the first step in the measurements. The measured values of mass attenuation 
coefficients of phantoms stimulating muscle, adipose tissue and bone were in good 
agreement with theoretically calculated values. For the twenty two energy channels 
used in this study, the ratio of the measured and theoretically derived values for the soft 

20 tissue varied from 0.982 to 1.004 and the average percent difference was 0.91%. For 
the adipose tissue, the ratio of the measured and theoretically derived values varied 
from 0.994 to 1.011 and the average percent difference of the measured values was 
0.54%. The ratio of the measured and theoretically derived values for bone varied from 
0.989 to 1.013 and the average percent difference of the measured values was 0.81%. 

25 Considering the tolerance in the thickness of the phantom measurements, 0.02 mm, all 
measured mass attenuation coefficients were within an acceptable limit. 

All the assembled phantom components simulated the values typically 
encountered in measurements of bone mineral in the spine and femoral neck. 
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For the measurement of the primary spectrum, the number of photons detected in 
different energy channels of the continuous x-ray spectrum varied from 110 s to 910 5 . 
The number of photons transmitted through different phantoms varied from 1-10 3 to 
4T0 3 . The calculated average percent difference between the calculated sum of the 
products of the mass attenuation coefficients and areal densities of the phantoms and 
measured fractional logarithmic transmission varies from 0.21% to 0.63%. There are 
several sources of error in the measurements of transmitted spectrum, but the most 
fundamental are the statistical fluctuation in transmitted intensity. 

When the photon absorptiometry method is used with less than three photon 
energies, the approximation of uniform adipose tissue distribution at the measurement 
site is used. The accuracy error, which results from the dual energy methodology, has 
been reported as high as 10%. (Michael, et al, "Monte Carlo modeling of an extended 
DXA technique." Phys. Med. Biol 43, 2583-2596, 1998). The reported precision of 
DXA is around 1%, but the results of bone mineral density obtained by different three 
manufacturers of dual x-ray absorptiometry equipment: Norland, Lunar and Hologic 
vary substantially. 

An alternative method of analyzing the three component body without the fat 
approximation is to use three different photon energies. As stated by several scientists, 
the triple energy method is less widely used, because the method is complex and 
requires a larger number of photons to be detected. {See Farrell et al. 9 1989, Kotzki et. 
al., 1991) This is the consequence of the small difference in the mass attenuation 
coefficients for muscle and adipose tissue in diagnostic energy range. The present 
study method used the continuous x-ray spectrum, measured using a multichannel 
analyzer and three component phantoms. The measured numbers of photons in the 
transmitted spectra were approximately the same as the number of transmitted photons 
used during the dual energy absorptiometry method. The x-ray unit used in this study 
allows achieving a very low tube current (10 fiA). This was important during the 
measurements of the primary spectrum, to reduce the pulse pileup. The tube current 
can be substantially higher during the phantom measurements. 
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Although the percent differences between the true sum of the products of areal 
densities and mass attenuation coefficients and measured fractional logarithmic 
transmissions were relatively small, solving the attenuation equation using Cramer's 
rule resulted in solutions that were not physically realistic. However, the sum of the 
solutions always gave the value close to the sum of the areal density of the phantom. 
For the five phantoms used in this study, changing the energy intervals from 1 to 10 
keV, this difference was always less than 5%. This fact was used as the base for the 
estimating routine in a computer program. 

In clinical studies, the spectrum transmitted through the soft tissue adjacent to the 
bone of the patient will be first measured. For the selected two energies the areal 
densities of the muscle and adipose tissue will be calculated using dual energy 
equations. In the next step the spectrum transmitted through the soft tissue and bone 
will be measured. The fractional logarithmic transmission will be calculated from the 
spectra over selected energy range. Using the Cramer's method the sum of areal 
densities will be calculated. The values of the areal density of the bone will be 
approximated by subtracting from the calculated sum of areal densities the sum of areal 
densities of muscle and adipose tissue calculated over the area of soft tissue only. The 
fraction of each component as the fraction of the sum of areal density will be 
calculated. Knowing that the inhomogeneity in adipose tissue distribution can 
introduce the error in bone mineral density grater than 10% in AP imaging and even 
greater in lateral bone mineral measurements the range of the fraction of each 
component will be calculated by adding and subtracting around 20% to the previously 
calculated areal densities. The best fractions of each component will be calculated 
using the computer program described previously in this project. 

As previously stated, to solve the linear equation using Cramer's rule the perfect 
values of fractional logarithmic transmissions should be used. This requires a large 
number of photons in the primary and transmitted spectra. Accumulating the large 
number of photons in the primary beam is not a problem, since this spectrum can be 
obtained earlier and stored in computer memory, but obtaining a large number of 
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photons in the transmitted beam will deliver a substantial dose to the patient. A 
computer program introduced in this study takes a different approach from the triple 
energy methods introduced by other scientists. The program calculates the areal 
density of the phantom components as the products of the best fraction of that 
5 component and the sum of areal densities of the phantom. The best fraction is selected 
from the estimating routine by calculating the minimum difference between the 
measured fractional logarithmic transmissions and evaluated values of sum of the 
products of mass attenuation coefficients, fraction and sum of the areal densities of the 
phantom. 

10 For the five phantoms described above, ranging the energy interval from 1 to 10 

keV, the average percent difference between the calculated areal densities of bone and 
the true values of that densities range from 0.42% for phantom 1 to 6.31% for 
phantom 2. No correlation was found between measured bone mineral density value 
and adipose tissue thickness. In the phantoms assembled in this study, bone mineral 

15 density was increased from 1.084 g/cm 2 to 2.168 g/cm 2 (phantom 4), and adipose tissue 
from 2.901g/cm 2 to 6.769 g/cm 2 (phantom 3). The computer program written in 
specific aim three underestimated the bone areal density value for four of the phantoms. 
The system will require calibration with more linear changes in areal densities of 
phantom components. 

20 The spectral information can be obtained by one of two methods. The standard 

method involves direct measurements using a multichannel pulse height analyzer. This 
requires expensive specialized equipment, which must be calibrated. The alternative 
method involves inferring spectral information from attenuation measurements. 
Deriving the spectrum from the transmission through the patient and increasing 

25 thickness of aluminum can be implemented in a clinical setting. 

A system and method for using x-ray spectra to determine areal densities of 
bone, fat and soft tissue in a localized area of a patient has been described. Although 
the above-described embodiment of the invention relates to a medical use, the invention 
has a number of non-medical uses, such as monitoring industrial materials for quality 
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control purposes, screening for hazardous materials such as explosives or weapons in 
baggage and packages for security purposes, and other uses in which it is desired to 
detect areal density of a composition in a multi-component mixture. It will be apparent 
to those skilled in the art that various modifications and variations can be made in the 
present invention without departing from the scope or spirit of the invention. Other 
embodiments of the invention will be apparent to those skilled in the art from 
consideration of the specification and practice of the invention disclosed herein. It is 
intended that the specification and examples be considered as exemplary only, with a 
true scope and spirit of the invention being indicated by the following claims. 
WHAT IS CLAIMED IS: 
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